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ABSTRACT 

Assessments of the viability of computational methods and the computer resource requirements for the 
prediction of tail buffeting are made. Issues involved in the use of Euler and Navier-Stokes equations 
in modeling vortex-dominated and buffet flows are discussed and the requirement for sufficient grid 
density to allow accurate, converged calculations is stressed. Areas in need of basic fluid dynamics 
research are highlighted: vorticity convection, vortex breakdown, dynamic turbulence modeling for 
free shear layers, unsteady flow separation for moderately swept, rounded leading-edge wings, and 
vortex flows about wings at high subsonic speeds. An estimate of the computer run time for a 
buffeting response calculation for a full span F-15 aircraft indicates that an improvement in computer 
and/or algorithm efficiency of three orders of magnitude is needed to enable routine use of such 
methods. Attention is also drawn to significant uncertainties in the estimates, in particular with regard 
to nonlinearities contained within the modeling and the question of the repeatability or randomness of 
buffeting response. 
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INTRODUCTION 


Current and future military fighters are called upon to perform high angle maneuvers at elevated 
loading conditions and are experiencing structural fatigue problems due to dynamic buffet loads on aft 
components of the empennage, in particular the vertical tail. In the past, steady maneuver loads on the 
aircraft have been predicted based upon a combination of model test data and analysis while dynamic 
buffet loads could only be estimated from model tests by means of suitable scaling procedures. For 
the most part, current military aircraft structural designs were accomplished prior to die realization of 
the extensive high angle maneuvering time the aircraft would accumulate and the resulting damage the 
structure would sustain. This has led to a desire for improved prediction methods for such dynamic 
buffet loads, which emerging computational fluid dynamics (CFD) methods might be expected to 
fulfill. The purpose of this paper is to provide an assessment of the issues involved in such CFD 
applications and the resources which would be needed to do such computations. 

Figure 1, from reference 1, illustrates the concepts of maneuver and dynamic or buffeting components 

of a representative structural response quantity, y B (t) . The maneuver component is the very low 
frequency, essentially steady state (D. C.), mean or bias of the signal. All loadings at frequencies 
below several Hertz are of this nature and include loads such as those due to short period 

maneuvering. The buffeting component, y B (t), comprises all of the signal content at higher 

frequences. A measure of the level of buffet response is the root mean square value of y B , a (t). 
The formula for a (t) given in the figure is an estimate of the square root of the area under a plot of the 

power spectral density of y B versus frequency. Fatigue and accumulated damage predictions based 

on assumed numbers of stress cycles are directly related to a and its accurate prediction will be critical. 
In order to calculate the buffeting structural response it is necessary to predict accurately the unsteady 
aerodynamic loading (buffet) causing the response. This in turn requires the accurate prediction of the 
source of the unsteadiness in the airflow. Details of buffet airflows on specific aircraft will be used to 
illustrate these issues. 

Configurations with Stalled Wing Flow, - Figure 2, from reference 12, sketches the buffet flow for the 
F-15 airplane model. Comments have been added to draw attention to the flow modeling issues. The 
F-15 wing does not develop intentional vortex lift augmentation and at 22 degrees angle of attack, 
where vertical tail buffet intensity is greatest, experiences essentially stalled flow over the outer wing 
panel with some underlying separated flow spilling over the inlet. The wing flow separates in the 
vicinity of the rounded leading edge in a highly unsteady manner. The subsequent flow convecting 
downstream past the vertical tails has an overall, average, diffuse vortical flow structure with higher 
frequency, unsteady swirling flow superimposed. Little is known about the dynamics within the 
separated region other than steady, average velocities and pressures. Low speed wind tunnel flow 
visualizations have been made. 

Configurations with Intentional Vortex Lift Enhancement. - Figure 3 shows details of the F-18 airplane 
which was designed to utilize the vortex system developed by the sharp leading edge extension for 
enhanced stability at high angles. Again, the comments draw attention to flow modeling issues which 
need to be addressed in any analysis of tail buffeting. The sharp edge forces the flow to separate at 
angles above a few degrees and the separated flow rolls up to form a vortex which is convected 
downstream by the mean airflow. For moderate angles, the vortex can remain intact for large distances 
downstream of the trailing-edge, although it does experience dissipation. For larger angles, the vortex 
can experience a fluid dynamic instability, termed breakdown or bursting. The bursting location 
moves upstream with increasing angle, reaching the trailing-edge at an angle a = 15-25 degrees. The 
burst location and the flow downstream is highly unsteady. Little is known about the dynamics within 
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the burst region other than steady, average velocities and pressures. A few measurements of rms 
values have been made. Water tunnel, low speed wind tunnel and flight flow visualizations have been 
made. 

The sections to follow will describe: experimental knowledge of vortex-dominated and buffet flows; 
the fluid dynamic equation levels available for calculation of such flows; estimates of the computational 
resources needed; and the current status of CFD applications. 


VORTEX-DOMINATED AND BUFFET FLOWS 

Experimental research in this area has been pursued separately by the aerodynamics and structures 
communities until recently. Reference 1 summarizes the traditional methods used within the structures 
community in correlating wind tunnel and flight buffeting response while reference 2 gives a recent 
assessment of dynamic loads due to flow separation. References 3, 4, 5 report buffet testing on the 
F/A-18, F-15, and F-l 1 1 TACT aircraft. Figure 4, from reference 3, gives buffet pressure frequency 
responses for the F/A-18 low speed wind tunnel model showing a marked effect of angle of attack 
upon the frequency for peak amplitudes. Reference 3 also shows this peaking to be a function of 
airspeed and a key to the tail buffet response via a "tuning" with tail vibration modes. This peaking in 
strake vortex flows is believed to be due to convection of patches of coherent swirling flow within the 
burst vortex. 

Reference 6 documents recent unsteady pressure and flow visualization tests on an oscillating straked- 
delta wing model similar in planform to the F-l 6. Figure 5, from reference 6, summarizes the main 
features of the two vortex systems (strake and wing vortices) for this model as evidenced by the lift 
versus angle of attack (AOA). Straked wing configurations like the F-l 8 and F-l 6 show the nonlinear 
vortex lift augmentation seen between 8-18 deg. and the downward break in the lift curve slope caused 
by the the onset of burst vortex flow over the planform. The maximum normal force coefficient occurs 

at a » 35 degrees. At higher angles the flow is fully separated and the value of C n falls. 

While the above research efforts focused on the flow unsteadiness in order to gain understanding of 
buffeting structural response, parallel efforts by the aerodynamics community have been directed at 
understanding the physics of separated vortical flows and the gathering of data bases for the validation 
of CFD codes. The experiments have involved detailed measurements on static, rigid models. Most 
of this basic research has studied vortex flows about idealized shapes; for instance, highly swept delta 
wings, sharp leading-edges, etc. Reference 7 summarizes results from the International Vortex Flow 
Experiment, a joint program studying vortex flow development on a 65 degree cropped delta wing. 
Tests were conducted between Mach numbers of 0.4 and 4.0. Both sharp and rounded leading-edges 
were tested for validation of Euler and Navier-Stokes codes, respectively. Also, sweep effects were 
studied with a 55 deg. swept wing and configuration effects studied with the addition of a canard. 
Reference 8 summarizes another extensive series of basic wind tunnel tests including an aspect ratio 
1.0, 75 deg swept sharp-edged delta wing and a double-delta (80 deg/60 deg) wing. Data bases from 
these tests include total forces and moments, surface pressures, surface flow visualization and flow 
field surveys via intrusive pressure probes (which disturb the flow being measured). 

Recognizing the need for higher quality flowfield data for code validation, recent tests 9 » 10 have 
utilized nonintrusive three component laser Doppler velocimetry (LDV). Reference 9 gives data for the 

75 degree delta wing at a = 20.5°. Reference 10 studied in detail the bursting of vortices generated by 
a 75 degree delta wing. This is also one of the first instances in which primary attention was given to 
unsteady features; root-mean-square (rms) velocity components are reported in addition to the mean 
velocities. Reference 10 also focuses attention to the boundary of the burst region just downstream of 
the burst location as the source of strong fluctuations. 
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The most recent wind tunnel investigations of vortex flow are beginning to show a merging of the 
interests of the two communities upon the problem of tail buffeting. Most of the earlier wind tunnel 
studies of vortex flows are not germane to this issue for a variety of reasons: i.) vortex systems were 
steady/stable, ii.) unsteadiness was not addressed or measured, iii.) idealized configurations (sharp- 
edges, highly swept wings, etc) were tested. Reference 1 1 describes LDV surveys of the flow over a 
YF-17 model in a low speed wind tunnel. Figure 6, from reference 11, shows mean and rms 

component velocities for a = 25 deg at a station just ahead of the vertical tail. The strake vortex has 
burst ahead of this station and the mean streamwise velocity component, u, shows a region of reversed 
flow centered on the vertical tail location. The maximum rms velocity fluctuations reach levels of 
approximately 40%, 35%, and 30% for the u, v and w components and mean crossflow angles at the 
tail vary from -30 degrees at the root to +15 degrees at the tip. The flow over the F-15 aircraft at high 
angles is different in not containing such high energy concentrated vortices (see fig. 2). Figure 7, 
from reference 12, shows mean crossflow velocity components at the tail station measured with LDV 

in a low speed wind tunnel and for a = 20 deg. Mean crossflow angles at the tail varied from +24 
deg. at the root to -17 deg. at the tip. Contrary to the burst vortex flow pattern shown in figure 6, no 
axial reversed flow is seen in the F-15 case. Also, reference 12 comments that hot-wire pressure data 
shows no high-frequency peaking like that shown in figure 4. However figure 19 of reference 4 
indicates the presence of dominant frequencies in the surface pressure spectra. 

This survey of wind tunnel testing of vortical and buffeting flows reveals the scarcity of data bases 
suitable for the validation of CFD calculations of buffeting flows. On the one hand, the aerodynamics 
community has been heavily occupied understanding and documenting capability for idealized low 
speed vortex flows. On the other hand, the buffeting and fatigue problems being experienced by 
current high performance fighter configurations are due to much more complicated flowfields. Vortical 
and separated flows from multiple aircraft components interact to produce the buffet environment. 
Also, compressibility and transonic effects at high subsonic speeds and high angles are largely 
unknown. Figure 8, from reference 13, is an example of the complex surface flow patterns that can be 
involved in realistic high performance aircraft buffeting. This surface oil flow for a high speed deep 
buffet penetration condition is a reminder that most of the experimental vortical flow research 
mentioned above has been for low speed conditions. The appendix to reference 7 discusses recent 
results on vortex bursting at high subsonic speeds, indicating that there is much that is not known 
about such flows. 


FLUID DYNAMIC MODELING 

The dynamics of fluid flows can be studied at many different levels and the complexity of the fluid 
flow modeling chosen must be carefully chosen to suit the application. Newsome and Kandil^ 
provide a concise summary of the governing equations and solution algorithms being applied to 
vortical flow aerodynamics and also provide a survey of physical modeling issues and computational 
results. The following discussion will emphasize the resources required for accurate computations. 
Thus, grid resolution for accurate solutions, the total number of grid points, the allowable time step 
size and computer speed will be the key players. 

Navier-Stokes Equations. The Navier-Stokes (NS) equations for a compressible fluid are derived from 
the conservation laws of mass, momemtum and energy. This yields five equations relating seven 

unknowns: density, p ; pressure, p; internal energy, e; temperature, T; and three velocity components 
u, v and w. The two additional relations needed to make the system solvable are provided by an 
equation of state, such as that for a perfect gas. An additional quantity introduced in the analysis is the 

coefficient of viscosity, |i. These NS equations describe the unsteady, turbulent motion of fluids. 
However, the small size of the turbulent eddy flows which they describe and the required grid point 
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density for accurate resolution preclude their use for flows about complex aircraft geometries. To get 
around this issue and bring engineering problems such as the aerodynamics of wings and bodies 
within reach, averaging procedures are used. The small scale, very high frequency eddy flows are 
eliminated from the modeling by a time-averaging process which produces the "Reynolds-averaged 
Navier-Stokes (RNS) equations." These RNS equations are capable of describing the steady and 
unsteady aerodynamic flows encountered in aeroelasticity and buffeting. 

In aerodynamic flows, the effects of viscosity are mainly confined to thin regions, termed boundary 
layers, adjacent to the wings and bodies. Figure 9 shows the growth of a typical boundary layer 

velocity profile over a flat plate. For attached flows, the boundary layer thickness, 8, is typically only 

a few percent of the body length and the velocity at the edge of the layer, U e , is that of the outer 

(external) inviscid flow. The viscous fluid nearest the surface must satisfy the no-slip boundary 
condition (U = 0). The viscous sheared flow within the boundary layer is responsible for the 
production of vorticity which in turn generates lift. The nondimensional Reynolds' number 


Re = p U 1 /p 


typifies the relative importance of viscosity on the flow. The larger the Reynolds number, the thinner 
the boundary layer, 5. Typical values of Re, based upon mean aerodynamic chordlength, are less than 
one million for low speed wind tunnel models and 10-50 million for full scale aircraft flight conditions. 

Real boundary layer flows over aerodynamic surfaces contain regions of laminar, transitional and 
turbulent boundary layers as illustrated in figure 10. Flow near a wing surface near the leading-edge 
stagnation point is nearly always laminar (purely sheared flow) but usually undergoes transition (i. e. 
the laminar flow experiences an instability) to turbulent flow in a short distance. The averaging 
performed by the RNS eqs. treats this addition of turbulence to a laminar flow by introducing an "eddy 

viscosity" coefficient, p t , giving the total viscosity as 

p = p. + p, 

“ “lam “l 


The introduction of eddy viscosity, p t , in the RNS eqs. requires that the effect of turbulence within 
the boundary layer profile be modeled. This process of turbulence modeling has been the source of 
extensive research (see reference 14). However, all popular turbulence models have been calibrated 
against equilibrium, wall-bounded shear flows (Fig. 9). At the present time there is no turbulence 
model that is applicable to a rotational, off-body flow (free shear flow), such as a vortex, either burst 
or unburst. 

The expense involved in RNS eq. calculations is due to the requirement for sufficient grid density to 
accurately resolve the boundary layer shear flow. This requires 10-20 grid points across the boundary 
layer. Computer resources are further impacted by the fine grid cells needed by the RNS eqs in the 
boundary layer since this affects the allowable time step size. To maintain numerical stability with 
such small grid cells, allowable time step sizes result in several thousands of time steps for each cycle 
of unsteady motion. This can be as much as an order of magnitude more time steps per cycle than 
required for the structural dynamic response calculation. 


Leading-Edge Vortex Flows. Figure 1 1 indicates the geometry of concentrated leading-edge vortices 
generated over a sharp-edged highly swept delta wing at a high angle of attack. The sharp leading- 
edge causes the flow to separate, creating a thin free shear layer (termed the feeding sheet) through 
which there are high gradients of the amplitude and direction of the fluid velocity. This viscous 
sheared region generates vorticity which is fed into the "primary" vortex core above the wing. These 
concentrated leading-edge vortices can contain high energies, with axial core velocities 2-3 times the 
freestream velocity. Reference 10 gives details of LDV surveys of such low speed vortices, showing 
that there is a.) an innermost, rotational, viscous subcore, b.) a middle region of inviscid, rotational 
flow surrounding the inner core and c.) an outer inviscid, irrotational flow. 

Figure 1 1 indicates a few of the many other details which may be present: reattachment of the flow 
inboard, followed by spanwise flow across the surface. The reattached spanwise flow on the surface 
creates its own boundary layer which can itself separate and form a small "secondary" vortex under the 
primary vortex. On the surface underneath each of these vortices, the high flow velocities induce 
suction pressures, adding to the lift generated by the wing. 

Leading-edge vortices do not persist in the flow indefinitely. Multiple vortex systems, like the two 
vortices in figure 11, mutually affect each other in downstream regions where they can intertwine with 
each other as they gradually dissipate due to viscous effects. Also, at higher AOA, the individual 
vortices experience instabilities or breakdowns. Table I, from reference 15, categorizes regions of 
vortical flow, terming this region of vortex instability Region III. Figure 12, from reference 16, 
shows two modes of instability occuring simultaneously on a wing; the "spiral" breakdown on the top 
and "bubble" breakdown on the bottom. Both forms are highly unsteady both in the location of the 
burst point and in the flow downstream. Both spiral and bubble vortex breakdown are characterized 
by an axial stagnation point and a limited region of reversed axial flow. The spiral form of breakdown 
occurs in three successive stages; a sudden deceleration of the fluid moving along the core, an abrupt 
kink of the core which then turns with a whirling motion that persists for a few turns and a breakdown 
to large scale turbulence. Particularly at flight Reynolds numbers, leading-edge vortex breakdown is 
an unsteady phenomenon characterized by the whirling motion of the spiral form or the oscillation of 
the bubble form about a mean position in addition to the downstream coherent unsteadiness. Figure 
13, from reference 37, gives a conceptual cross sectional schematic of the spiral bursting. The 
individual coils maintain a large scale coherency for a limited distance and, convected by the freestream 
flow, will give rise to a characteristic frequency content as observed at a fixed location on the aircraft 
(e.g. the vertical tail). Such frequency peaks have been noted in buffet tests in the F/A-18 (reference 
3) and in low speed wind tunnel tests Figure 14, from reference 10, shows the region of strong 
rms velocity fluctuations just downstream of the bursting location and forming a shear layer conical 
boundary of the region of reversed axial flow. Numerical analysis^ of high speed film images 

established that the convection velocity of the spiral cores was close to U„/2 and that their frequency 
was close to 130 Hz., corresponding to a Strouhal number St = fd/U = 0.7 (d is the maximum 
diameter of the recirculating bubble). 

Reference 15 states "there is no clear understanding of the mechanisms leading to vortex breakdown. 
Essential elements seem to be a region of high total pressure loss within the vortex core and an 
externally imposed axial pressure gradient.. ..despite numerous attempts, a generally accepted 
theoretical description of vortex breakdown does not exist." Alternative explanations of breakdown 
mechanisms fall into three groups: hydrodynamic instability theories, stagnation point theories and 
wave-motion theories. Note also that the properties describing these vortex flows prior to bursting are 
frequently regarded as having smooth variations. For instance, the free shear layer and rotational 
vortex cores are regarded as smooth thin sheets and regions. However, upon closer observation more 
detailed dynamic structures are seen. Figure 15, from reference 17, illustrates intricate details seen in 
laser light sheet profiles taken before and after the breakdown. The condition is for an 85 degree delta 

wing at a = 45 degrees. The edge of the core and the shear layer exhibit localized patches of vorticity 
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forming giving a "lumpy" texture to the flow. What level of detailed simulation of such features 
required for accurate buffeting calculations is largely unknown. 


is 


This survey of leading-edge vortex flows shows that much is known about such flows for highly 
swept wings ( A> 60 deg.) and, in particular, sharp-edged wings. Much less is known for 
conditions typical of main wing panels of current fighter aircraft (i.e. A = 40-50 deg. with rounded 
leading edges). In these cases, it is known that there is no clean pattern of flow developing with 
increasing angle of attack. Thus, the flow contains features of both leading-edge vortex flow 
formation and stalled flow. Vortex flow separations may develop over partial regions of the span and 
be intermittent both in time and spanwise location. Vorticity generation for such conditions will not be 
easily treated as confined to a single shear layer emanating from a fixed location on the wing as in 

fianr#* 11 ® 


pli l.gr Equations, If the terms involving the viscous shearing stresses are deleted from the RNS eqs. the 
inviscid Euler equations result. This model does not require very fine grids near bodies since the 
viscous boundary layer is not resolved. The no-slip boundary condition is replaced by a condition that 
the flow direction at the surface be parallel to the surface. These equations are capable of modeling 
flows with strong shocks and flows containing vorticity (rotational flows). Several issues related to 
the use of the Euler equations for vortical and buffeting flows which should be understood involve the 
sources of vorticity production and instability mechanisms in free shear layers and concentrated 
vortices. The former issue is probably most significant in the "forward” portion of the flow-the 
orebody and wing leading-edge regions. The latter issue predominates in the downstream region 
where vortex bursting and large scale swirling motion are encountered. There have been numerous 
publications showing that the CFD solutions of the Euler equations for sharp leading-edge delta wings 
do produce vortex flow solutions (for medium sized grids) which are in surprising agreement with 
important features measured experimentally. Also, medium grid solutions have been shown to predict 
vortex core axial stagnation points and regions of reveresed flow, thus seeming to simulate vortex 
bursting. The controversy oyer the accuracy and convergence of such solutions for sharp and round 
leading-edge wings is described in detail in reference 15. Reference 18 is representative of Euler 
equation calculations for such flows while reference 19 concludes that, for blunt/rounded leading-edge 
wings, Euler equations methods do not give satisfactory results. 


STRUCTURAL DYNAMICS SIMULATION FOR BUFFETING RESPONSE 

Traditional structural loads analysis methods distinquish between maneuver and dynamic or buffet 
loads whereas the CFD methods under discussion are, in principle, able to compute the total load 
directly. This does not come without a price however, since the very low frequency (maneuver) loads 
must be computed with the same restrictive time step needed for the highest frequency components. 
Since these low frequency loads can be treated by other methods, a low frequency cutoff limit of 5-10 
Hz. will be assumed for buffeting calculations. Also, buffet load spectra typically involve the lower 
/^ru- UenCy V1 j llon mo ^ es so titat an upper cutoff frequency may be assumed in the range 40-100 hz 
(Ihis paper does not address acoustic resonances/fatigue which occur at higher frequencies and 
involve local structure panel response.) 


wACTD^Kn 0 °j l , he f ai1 structure will be assumed given by a finite element model (e.g. 

NASTRAN). The model will be truncated to the 5-20 lowest frequency modes with each mode 
simulated by an ordinary differential equation of motion. Buffeting response calculations mav be 
performed in several ways. The buffet load environment at the tail is assumed to be calculated by the 
LrD code as time histories of integrated loads (lift, moment, generalized aerodynamic forces (GAFs)). 
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- One- way-coupling. Here the modal equations are integrated in time along with the CFD code giving 
the buffeting structural response and loads. However, it is assumed that the tail motions do not affect 
the aerodynamic flow environment at the tail. Thus, the buffet loading on the tail is treated as an 
external forcing function. Since the truncated mode structural equations are linear, frequency response 
analysis can be used to calculate buffeting response quantities such as power spectral densities and rms 
values. It is imperative in this approach to account, in some fasion, for aerodynamic damping. Total 
damping is typically low and the major portion of the buffeting response is concentrated in narrow 
frequency bands centered on the structural mode resonances. Aerodynamic damping 1 is proportional 
to dynamic pressure and usually dominates the structural damping. Integrated buffeting response 
quantities such as rms bending moments are strongly dependent upon the aerodynamic damping. This 
one-way-coupling method is similar to the Rigid Tad Pressure Method of reference 3. 

- Two- way-coupling. In this most complete simulation, the buffet load environment drives the 
structural equations of motion and the resulting structural motion influences the CFD flow calculation. 
At each time step, the computational grid about the tail must be moved to maintain alignment with the 
vibrating tail surface and the updated surface location is used as the new boundary condition by the 
CFD code. 

The two-way coupling method has been used for aeroelastic calculations with CFD codes for attached 
flow conditions. 3 ® However, neither of these methods has been demonstrated for buffeting response 
calculations using CFD codes. 


COMPUTATIONAL RESOURCE REQUIREMENTS 

We turn now to the estimate of the computer resources which will be required for vertical tail buffeting 
calculations. The goal will be to obtain at least a ballpark estimate of the computer memory size and 
job run times. Until recently, the former item was the pacing item in CFD calculations, with the 
allowable number of grid points being restricted by computer memory sizes. The currently available 
supercomputers, with memories ranging from 32-256 million words of core memory, have made job 
run times the pacing item as the estimates below will show. 

The CPU run time for a CFD calculation can be estimated from the relation 


where 


0) 


N s t = number of time steps in the calculation 

N gp = number of grid points 
x = computational time per grid point per time step 

Algorithm Speed -The parameter, x, is a common measure of the speed of an algorithm and typically 
ranges from 10 - 100 microseconds per grid point per time step. This assumes cun , ent computer 
operational processing speed capability on the order of 250 million floating point operations per second 

(MFLOPS, typical of the CRAY 2). The value, x = 10, is typical of "explicit" algorithms which are 

commonly used for steady-state CFD calculations. Implicit algorithms obtain values of x = 20 - 100 
since they require more complex programming. However, they are favored for time-accurate 

calculations since they allow larger time steps.. A value of x = 40 will be assumed herein. 
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Number gf Computational Steps -The parameter N f t depends upon the time step size, At, and the 
total real time length necessary for the calculation, T t ot (to be distinguished from total CPU run time, 
T ) 

cp U /• 


N st=T lol Mt 


( 2 ) 


The maximum time step size, At, is limited by the numerical stability of the algorithm and by the 
required accuracy for the aerodynamic results. To perform a buffet calculation, the frequency 

bandwidth of interest must be specified. Let f min and f m>x designate these frequency limits. 
Corresponding maximum and minimum periods of oscillation are then given by 


T . = 1/f 

min max 


T = 1/f 

max 1M min 


For the F-15, let f m . n - 10 Hz. and f max =40 Hz. Now, the time step. At, is chosen as 


At = T . /N , 

min s t /cy 


(3) 


where N ft/cy is the number of computational steps per cycle of oscillation. For Euler codes this 

value is set by the required accuracy for the unsteady calculations and is in the range 100-200. For 
Navier-Stokes codes, smaller time steps are required due to algorithm stability and values of 1000 and 
higher are needed (due to the much smaller minimum grid size used in order to resolve the boundary 

layer.). In the examples to follow, N $t/cy values of 200 for Euler calculations and 1000 for Navier- 
Stokes calculations will be used. 

The total real time length, T 1 0 , , required for buffeting calculations will be set by the number of cycles 
of oscillation, N c y , at the lowest frequency, f m ; n , in order to obtain converged results. 


T = N * T 

lot ^cy 'max 


(4) 


Since the buffeting flowfield is inherently nonsteady and nonperiodic, convergence must be measured 
in some statistical sense such as the attainment of stabilized power spectral densities of selected 

structural response parameters. A reasonable range of N c is 10-50. A value of N c =10 will be 
assumed herein. cy 

Summarizing, the total number of steps is given by 


( 5 ) 


N s t = T tot Mt 

= (N cy *T max )/(T min /N st/cy ) 
= (N cy /f ra . x )/(l/N st/cy *f m#x ) 

= N cy*N gt/cy *(f max /f inin ) 


For the values selected for the case of the F- 15 this gives 

N st = (10) * (200) * (40/10) = 8,000 steps (Euler) 

= (10) * (1000) * (40/10) = 40,000 steps (Navier-Stokes) 

Number of Grid Points - In this section, a first approximation will be developed for the number of grid 
points needed for a buffeting calculation of the F-15. First, the current CFD experience with solution 
convergence for isolated wings will be noted. Then attention must be given to the different nature of 
buffet airflows. These two features will be referred to as "local accuracy" and "remote accuracy." 
Figure 16, from reference 20, shows the near field of a typical grid used for CFD calculations about an 
isolated wing. Only the left half-wing is shown since symmetry is usually assumed to cut the total 
number of grid points required in half. The grid point density is high near the surface of the wing 
(dark area) where high accuracy is desired for steady calculations, and the density is decreased rapidly 
away from the wing where lower flow gradients are anticipated and where less accuracy is acceptable. 
Calculated pressures on the wing surface are used to compute the wing's lift, moment and drag. Thus, 
the accuracy of these "local" flow quantities is of most interest. On the other hand, for buffet 
calculations, the accuracy of the flow quantities (e.g. velocities and pressure) at the "remote" tail 
location will be of paramount interest. For local accuracy with an Euler code on a half-span wing, 
reference 21 indicates that the following number of grid points results in spatially converged local flow 
results: 


193 grid points in the stream wise direction 
33 grid points in the direction normal to the wing 
41 grid points in the span wise direction 

The total number of grid points is N gp = 261,000. Reference 20 gives similar information for a 
Navier-Stokes code calculation: 

193 grid points in the streamwise direction 
65 grid points in the direction normal to the wing 
49 grid points in the spanwise direction 

The total number of grid points is N gp = 414,000. Note that the only coordinate direction with 
significantly higher grid density is that normal to the wing. This is due to the use of the "Thin-Layer 
Navier-Stokes" equations (TLNS) wherein it is recognized that the viscous shearing forces in the 
direction normal to the wing predominate in boundary layer flow. Accordingly, the viscous terms in 
the streamwise and spanwise directions are not treated and no additional grid point density in these 
directions is called for, leading to a very large reduction in the number of grid points for a viscous 
calculation. Comparing the two grids, it can be seen that approximately 20 additional grid points in the 
normal direction are needed to resolve the viscous boundary layer. Thus, the TLNS assumption 
allows RNS-like results with only a 30-40% increase in cost over Euler codes. We note that there are 


1 1 



approximately 100 streamwise grid points on the wing itself, leading to 100 X 20 X 49 = 98,000 grid 
points in the boundary layer of the wing. If the use of the thin layer approximation cannot be justified 
resolving viscous regions becomes prohibitively expensive. Thus, doubling the resolution for general 
sheared flows using the RNS equations requires an eightfold increase in the number of grid points 
(due to increased grid sizes in all three coordinate direction, 2^). 

Turning now to buffet computations, the issues of unsteadiness and the convection of vorticity must be 
addressed. Now the grid stretching shown in figure 16 cannot be used in the vicinity of the wing- 
fuselage-tail since remote accuracy is required (the effect of disturbances originating in the wing region 
must be accurately resolved at the tail). For this purpose, the computational domain will be treated as 
three regions. In Region I, the boundary layer region of the wing, fuselage and tail, TLNS grid 
densities will be used to resolve the boundary layers. In Region II, the near-field of the aircraft 
containing unsteady separated vortical buffeting flow, grid densities adequate to resolve and convect 
the flow for the frequency bandwidth of interest will be assumed. Region HI comprises the outer flow 
field (inviscid and irrotational) necessary to capture the correct global flow. Thus the total number of 
grid points is 


XT XT 1 .XT 11 , xtIH 

N = N + N + N 
8P gp gp gp 


(7) 


The number of grid points required for Region III will be a small fraction of the total number of grid 

points and will be assumed to be = 50,000. The number of grid points for the boundary layer 
about each body component will be taken from the TLNS boundary layer grid density cited above:^0 


Ng p = (100 X 20 X 49) X 2 X N bod 


= 196,000 X N, . 

bod 


where N bod = the number of component bodies and the factor of two is used for full span modeling. 
For a wing-fuselage-tail configuration, N bod = 3 and N' gp =588,000. 

Figure 1 7 indicates the near- field Region II for the F- 1 5 aircraft where grid density must be adequate to 
accurately resolve and convect the unsteady features of the flow for the bandwidth of interest Since 
most of the energy involved in separated flow buffet is contained in the frequency range of the lower 
smictrual vibration modes, roughly 10-100 Hz., CFD methods must accurately predict the dynamics 
of the fluid at these frequencies wherever they may occur in the flow. In this section, it will be 
assumed that the correct amount of unsteady vorticity has already been injected into the flow (e.g. via 
the viscous boundary layers, Region I). Hence, in Region II the Euler equations will be assumed 
adequate to accurately convect the vorticity. The following assumptions are made in order to estimate 


the near-field Region II shown in figure 17 should have a grid density adequate to resolve buffet 
frequency components up to f m a x which will be assumed to be 40 hz. 



fluid disturbances will be assumed to propagate in the stream wise direction at the freestream 
speed, U „ . This leads to a minimum spatial wavelength in the x-direction of 


W„= U -' f n..x ft/cycle (8) 

- to accurately resolve streamwise fluctuations with wavelengths of X m . n _ x , N gp/w j grid 

points per wavelength will be assumed. Typical values for N gp/wl are 50-100 ( N gp/wl =50 
will be assumed herein.) 


Ax = 


X . /N . . 

min-x gp/w 1 


(9) 


= MN„ /wl Xf mn ) 

- for the large flow angularities found in buffet flows, it will be assumed that the spatial 
resolution needed for the crossflow directions will be similar to that for the streamwise resolution. 

Ax = Ay = Az (10) 

- assume a freestream Mach number of M = 0.5, giving U„, = 500 fi/sec. 

- for the dimensions shown in figure 17, these assumptions lead to the following grid parameters 

X . = k min _ v = X . = U„/f =500/40 =12.5 ft/cycle 

min-x min-y min-z 00 max J 

Ax = Ay = Az = X / N . . = 12.5/50=0.25ft 

J min-* gp/w 1 

and the number of grid points in the near-field Region II shown in figure 17 is calculated to be 
approximately 1.5 X 10^ 

Thus from eq. (7) the total number of grid points for an F-15 buffet calculation is 

N gp = 588,000 + 1,500,000 + 50,000 
= 2,138,000 


Finally, the total CPU run time for one F-15 buffet calculation (full-span), assuming that the TLNS 
equations are used in Region I and that this determines N st , is estimated to be 


T cp u = 40,000 X 2,138,000 X (40X1 O' 6 ) / (3600 sec/hr) 
= 950 hrs 


If an Euler equation calculation is assumed ( N 


I 

gp 


= 0, no viscous boundary layer, N s t 


= 8,000) then 


T c P u = 8,000 X 1,550,000 X (40X10-6) / (3600 sec/hr) 
= 138 hrs 


Equations 1-10 are summarized in Table II. While these extremely long run times would seem quite 
impractical, they are likely to indicate the range within which solutions with sufficient accuracy to be 
useful for design decisions will be found. On the one hand, the shorter run time obtained assuming 
the Euler equations represents an absolute minimum amount of modeling, evidenced by the lack of any 
physical model for vorticity generation. At the other extreme, the excessive run time of the TLNS 
example should be capable of being reduced by a number of devices such as; restricting the extent of 
boundary layer modeling to only those regions essential to simulating unsteady flow separation 
(wing/strake leading edges, etc), implementing time-accurate local time stepping, and other 
possibilities. 


CURRENT CFD APPLICATIONS AND CAPABILITIES 

The above preliminary estimates of computer resource requirements for buffet calculations will now 
be juxtaposed with the current capabilities of CFD. Representative published applications of CFD 
methods will be discussed for: complete aircraft configurations near design conditions; vorticity 
convection; and high angle, steady/unsteady, vortex-dominated cases. 

Complete Aircraft Configurations. It is not surprising that the most detailed aircraft geometry 
modeling has been applied for cases at low angles of attack, near design conditions, where the flow is 
attached or mildly separated. These conditions are most important in design and the codes can be 
expected to perform at their best due to good flow quality (steady, attached, thin shear-layers, etc). 
Thus, indications of the accuracy of results here will help to assess the readiness of the codes for the 
more demanding buffet conditions. Table III summarizes a number CFD applications for such cases. 
All of these studies used the TLNS equations and all implemented simple turbulence models. All used 
convergence acceleration devices (local time-stepping, multigrid, etc) which yield accurate results only 
for converged, steady flows. Except for reference 23, all made use of assumed symmetry in the flow 
to reduce by half the grid size (only one-half aircraft model is analyzed). Finally, the codes used in 
these studies were generally second-order accurate in space and first-order accurate in time. 

References 22-24 involve complete aircraft models at transonic speeds while reference 25 involves a 
detailed forebody/strake model at low speed and high angle. The latter is included due to its complex 
surface modeling. Fujii and Obayashi 22 modeled the W-18 transport configuration as a wing- 
fuselage and made calculations for three angles. The overall surface pressures compared fairly well 
with experiment except in the outboard wing region where elastic deformations were not accounted 
for. Their numerical algorithm allowed the calculations with 700K grid points to be obtained in 5-6 
hours. Flores and Chaderjian’s 23 study of the F-16A aircraft is one of the most ambitious applications 
to date. They modeled the wing, fuselage, tail, inlet, inlet-diverter, and the exhaust nozzle using 27 
grid zones for the half-airplane. The flow-through inlet was modeled including power effects. 
Comparisons with experimental pressures indicated that the wing leading-edge expansion was not 
adequately resolved and the wing shock location was off by 12% chord. The differences are suspected 
to be due to insufficient streamwise grid resolution. Doubling the grid size to one million points 
allowed the calculation of the full-span aircraft at five degrees sideslip angle. This also doubled the 



runtime to 50 hours. With regard to buffet calculations, it is interesting to note that the vertical tail tip 
vortex for this condition dissipated within one tip chordlength due to numerical dissipation (due to grid 
stretching downstream of the tip). Huband et al.24 studied the same F-16A (the inlet is faired-over) 

for a low supersonic Mach number. Their fine grid solution (N gp = 1,24 IK) occupied 59 million 
words of memory and required 40 hours of runtime. They obtained favorable agreement with 
experimental surface pressures but the wing leading-edge suction peaks were not correctly predicted 
due to lack of numerical resolution. 

Thomas et al.^5 studied the vortex flow patterns over the F/A-18 forebody and leading-edge extension 
(LEX) at low speed and for a = 30 deg. Both laminar and turbulent calculations were made and 
compared with low Reynolds number, low speed wind tunnel tests and with high Reynolds number 
flight tests. In terms of surface flow patterns, primary and secondary forebody surface separation 
lines were well predicted. The laminar calculations (Re = 740K) agreed well with the wind tunnel (Re 
= 200K) results while the turbulent calculations agreed well with the flight tests (Re = 10 7 ). The two 
calculations show significant differences in surface flow patterns indicating the effects of laminar 
versus turbulent flows. The downstream convection of these vortical flows under buffeting conditions 
can be expected to lead to transition and turbulence effects on buffet. 

In summary, ambitious applications of CFD to complicated aircraft geometries are being performed. 
Available studies are very encouraging in the overall agreement with experiment. Important areas are 
also being highlighted where additional grid resolution is needed to achieve local accuracy in such 
important features as suction peaks and shock locations. Computer runtimes for these cases, where 
remote accuracy is not required, are in the range of 5-40 hours (half-airplane). It is probable that 
accurate buffet calculations will require capability such as this for forebody and wing flows in order to 
generate accurate "starting" conditions for the convecting vortex flows. It is also probable that similar 
capability will be required to calculate the buffeting response of the tail to the oncoming buffet flow. 
Indeed, this is likely to be a more difficult problem than the calculation of the "start ing" wing flow due 
to the turbulent nature of the local flow at die tail (see figure 6). 

Vorticitv Convection . Any finite difference numerical algorithm must contain a dissipative mechanism 
for stability. Standard methods involve either the explicit or implicit addition of nonphysical damping 
terms proportional to a combination of second and fourth spatial differences of the unknown variables. 
The resulting CFD algorithms have been typically second-order accurate in spatial variation. 

For calculations of vortex flows, where the vortex is calculated as a part of the solution and not 
imposed on the solution, a significant problem is that conventional differencing schemes possess 
enough artificial dissipation to smear and dissipate the vortex very rapidly. In general, the coarser the 
finite-difference grid that is used to perform the calculation, the greater the distortion of the vortex. 
Rai^ti indicates that conventional spatially second-order accurate finite-difference schemes are too 
dissipative for calculations involving vortices that travel large distances. Studying the model problem 
of a two-dimensional, crossflow vortex being convected by the ffeestream, the effect of higher order 
spatial- and time-accurate differencing was evaluated. A commonly used method (Beam-Warming 
central spatial differencing and first order accuracy in time) was shown to dissipate more than 20% of 
its core pressure in only 5 vortex core radii of travel, making it not suitable for such calculations. 
Increasing the time-accuracy to second order is very effective, decreasing the pressure dissipation to 
20% in 45 core radii. Increasing spatial accuracy is also required to further decrease the numerical 
dissipation: a fifth-order upwind-biased spatial scheme with second-order time-accuracy produced only 
3% pressure loss in 45 core radii. To achieve these accuracies, there is also a lower limit on the 
number of grid points within the vortex core. Reference 26 indicates that an approximate value for this 
grid density is at least 8-10 grid points per core diameter of the vortex. 



No similar study of accuracy and grid density requirements for strcamwise vortices such a. ww 
edge vortices, has been made. Considering that the initial scale 
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and angle-of-attack upon the vortex flows for angles below bursting. A clear well-organized vortex is 
indicated for the 60 deg wing which moves inboard with increasing angle. The double-delta wing adds 
an additional vortex formed at the strake-wing juncture which interacts with the wing vortex, intertwining 
with it over the wing. With the cropped double-delta wing, an additional vortex appears due to the flow 

separating from the side edge and at a = 20 deg. all three vortices tangle together in the wing-tip region. 

Fujii and Schiff^O present calculations for the same 80-60 deg. double-delta wing using a time-accurate, 
compressible TLNS code. Results for angles of 12-35 deg. are shown. The lift coefficient 

agrees well with experiment up to a = 27 deg. where breakdown occurs. For the 850,000 point grid 

used the strake vortex at a = 12 degrees tended to loose strength rapidly because of numerical 
dissipation. For 30 and 35 degrees, unsteady vortex breakdowns were calculated. At 30 deg. a bubble- 
type breakdown was seen at x/c = 0.85, while a spiral-type breakdown occurred for 35 deg. at x/c = 0.6. 
Figure 21 of reference 30 shows the time history of the lift coefficient, for the latter case, containing an 
oscillatory component, presumably due to the spiraling bursted flow, with frequency f = 0.3 U / c. This 

gives a Strouhal number, based upon maximum burst region diameter, d, of St = 0.1. These calculations 
were repeated with the viscous terms deleted in order to simulate the Euler equations and the results 

differed significantly from the viscous case. For a = 12 deg. the two vortices didn't merge over the 

surface; for a = 30 deg. breakdown was not observed; while the a = 35 deg. case indicated breakdown. 
Reference 27 concludes that still better grid resolution is required to obtain quantitative results for such 
flows. 

Thomas et al.31 studied the 75 deg swept delta wing using a local time-stepping, compressible TLNS 
code with an 850,000 grid point mesh. Results for angles of 0-40 deg. are shown. The predicted 
maximum lift coefficient of 1. 10 at 35 deg. agrees closely with the measured maximum lift of 1.06 at 33 
deg. At 40 deg. a steady, bubble-type reversed flow region (due to breakdown) extending from x/c = 0.6 
to just downstream of the trailing-edge was observed. The experimental burst location for this condition 

is x/c = 0.4. Detailed comparisons with experiment for a = 20.5 deg show good agreement of minimum 
pressure coefficients under the primary and secondary vortices. Additional calculations for this wing are 

reported by Taylor et al.34. Comparisons with 3-D LDV flowfield measurements for a = 20.5 deg. 
show that the CFD code predicts the physics of the flow well, such as the position of the primary and 
secondary vortex cores and the position of the secondary separation line. However, differences in the 
magnitudes of pressure, velocity and vorticity in the region of the primary vortex core were noted. Krist 
et al.32 s how that these deficiencies are due to grid resolution in that region. Using a grid embedding 
technique to achieve high grid density there without the necessity of global grid enrichment, it is shown 
that the deviation from die experimental core velocity approaches zero as the minimum grid spacing in the 
core region is decreased. Accurate prediction of such vortex core quantities will be needed for accurate 
buffet predictions. 

The final entry in Table IV is that of Vadjak and Schuster^ who made calculations for a generic fighter 
configuration consisting of a sharp-edged strake/wing/fuselage. Low speed wind tunnel LDV flowfield 

data was available and crossflow velocity comparisons for a = 21 deg. appear to be good. A bubble- 
type reverse flow region, indicating breakdown, is observed at approximately x/c = 0.9. These 
comparisons for this high angle, vortex flow case are very noteworthy, particularly due to the complex 
geometry tested. However, the computational results, obtained with a grid of 101,000 points, cannot be 
regarded as converged. 

In summary, calculations of vortex-dominated flows about simplified highly swept wing geometries at 
low speeds are available. Early Euler equation results gave encouragement to their use for such flows, but 
recent detailed studies point to issues regarding their accuracy and convergence for realistic geometries. 



The thin-layer Navier-Stokes equations, both laminar and turbulent, are being used. The results capture 
key features of the flowfield such as primary and secondary vortices and surface pressure details 
Indications that vortex breakdown is being simulated have been published. However, all cases surveyed 
indicated the need for further grid refinement in order to achieve quantitative agreement with experiment. 
Current grid sizes of 250,000 to 10^ grid points are being used for half airplane modeling, leading to 
CPU runtimes of 2-25 hours. These grids are typical of those developed for attached shear flows with 
the highest gnd density near the body. This leads to inadequate density in off-the-body regions where 
concentrated vortices are located and mesh enrichment methods are being developed to address this 
problem. A fundamental problem is the lack of any turbulence model designed for dynamic free shear 
layers. 


DISCUSSION 


Equation (1) for the total CPU runtime for a buffet calculation can be rearranged to give 
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( 11 ) 


This expression indicates a strong dependence of T cp u upon f ma x and U M due to the grid density 
in the near-field. Region II. For this example, increasing f max from 40-60 Hz increases T cpu by a 
factor of 3.5 while halving the freestream speed, U M , from M = 0.5 to 0.25 results in a six-fold 

increase in T cpu . The F-15 buffet example indicates that approximately 1000 hours will be required 

for a TLNS code buffet calculation while a calculation with an Euler code would require 138 hours. 
The Euler result almost certainly does not have the necessary basic physical modeling and the TLNS 
runtime is unacceptable for analysis or design applications. In this section, the realism of the TLNS 
runtime will be discussed, approaches to reducing it will be mentioned, and needed comparisons with 
experimental buffet studies will be outlined. 


TLNS Runtime. Equation 1 1 assumes complete aircraft modeling and thus makes no assumption of 
symmetry. This is appropriate since asymmetries in the vortical flows over the wing panels should be 
anticipated (see Table I). Indeed, it should probably be assumed that the two vortex systems will 
interact dynamically with each other. The doubling of the runtime due to the doubled grid size must be 
accepted as necessary. 

Equation 11 also assumes that the global time-step is set by the most restrictive criteria - allowable 
time-step for stability on the smallest grid cell. This restriction is necessary in order to maintain time- 
accuracy with available codes. 


Equation 1 1 does not include any allowance for grid points to model free viscous shear layers away 
from the aircraft component surfaces. Regions that might be anticipated to require such extra gridding 
are: i.) concentrated vortex cores such as those shown in figure 3; ii.) free shear layer surfaces such as 
those shown in figures 11 and 15 and in the flow separations from rounded leading-edges (figure 2)- 
and iii.) conical surfaces of sheared flow on the boundary of the burst region just downstream of 
vortex burst locations, see figures 3 and 14. If these regions were steady, straightforward estimation 
procedures (e.g. 10 grid points across concentrated vortex core diameters, 10-20 grid points in the 
normal direction through thin free shear layers) seem to indicate modest increases in the number of grid 
points for a reasonable number of such viscous regions. Unfortunately, these regions are not steady 



and so a local grid enrichment strategy will not be successful unless provision is made for moving the 
dense grid region to track the evolving flow. Also, these regions will not be aligned with the 
coordinate system in such a way as to take advantage of a "thin-layer" assumption. 

If it is assumed that the oncoming freestream flow is steady, there must be a region where unsteadiness 
first occurs. We can anticipate that the unsteadiness will be due to local instabilities in the flow which 
are likely to be highly sensitive to a number of parameters. It is also likely that nonlinear mechanisms 
will be involved, which can lead to extreme sensitivity to initial conditions. This can call for very 
extensive calculations in order to cover reasonable ranges of conditions or to achieve statistically stable 
calculations (e.g. stable power spectral density estimates). There is also the possibility of chaotic 
response, in which there is no repeatability in the response. In this light, the assumption that the 

number of cycles of oscillation at f . be in the range 10 < N cy <50 might be very optimistic. 

The CPU runtime estimate for the TLNS equations needs to be reduced by approximately three orders 
of magnitude (from 1000 hours to approximately 1 hour) before such buffet calculations will be viable. 
Then the number of cases, the job turn-around time and the computational expense would be at levels 
that could support efforts leading to reliable engineering tools. These runtime improvements can 
probably be anticipated due to expected increases in computer speed and memory and increases in 
algorithm speed and efficiency. 

Efficiency Enhancements . There are several direct approaches which should be of use in reducing the 
CPU runtime estimates. Embedded grids and solution adaptive gridding strategies have already been 
mentioned wherein the evolving flow solution is searched for regions with high gradients. Additional 
grid points can then be added selectively in these regions resulting in good accuracy with many fewer 
total grid points. Another opportunity for reducing the number of grid points is the use of 
"unstructured grids." Structured grids, like that of figure 16, are constructed along orthogonal 
coordinate directions, resulting in quadrilateral grid cells. Adding grid points locally to such meshes 
results in additional grid points being added throughout the mesh. Unstructured grids are typically 
based upon tetrahedral cell elements and are not arranged along coordinate axes. By giving up the 
connectivity associated with arranging grid points along coordinate axes, the unstructured approach 
allows for high densities of tetrahedral cells (whose size can be arbitrarily selected) to model surface 
geometry in great detail and without requiring global grid enrichment. References 35 and 36 describe 
the development of these unstructured grid and adaptive grid refinement methods for aeroelastic 
applications. 

The use of multiple CPU processors working on the same problem, termed parallel processing, also 
has the potential of reducing the computational runtime. Reference 33 reports a speed up of a factor of 
10 with the use of vectorization and two parallel processors. Research is also being reported on time- 
accurate local time-stepping which may yield time savings by circumventing the requirement of using 
the most stringent allowable time-step size as the global time-step. 

While this paper has focused upon finite difference solutions of the flow equations, there are 
alternative solution methods which might have advantages for buffet calculations. However, only 
finite difference methods (and the closely related finite volume methods) have advanced to the stage 
where direct buffet computations appear feasible. Attention should always be given to tradeoffs 
between highly accurate, expensive solutions and inexpensive, semi-empirical solutions. In this vein, 
it might be anticipated that somewhat general guidelines will emerge in the course of the computational 
buffeting studies which will, without doubt, occur. 

Experimental Vortex-Dominated and Buffet Studies . This paper has discussed a number of areas in 
which measurements of vortical and buffeting flowfields are needed in order to validate or verify 
computational buffeting predictions. 


- Unsteady surface pressure tune-coirelations. Knowledge of the instantaneous surface pressure 
distribution over the vertical tail is required in order to calculate the tail buffeting response. Time 
correlated measurements of these pressures are needed along with off-the-surface flowfield 
measurements in order to begin to understand the interaction between the tail structure and the 
oncoming buffet flow. This constitutes extremely large data sets and presses the current limits of 
nomntrusive measurement techniques. 

Basic research into dynamic free shear layer turbulence models is probably vital to achieving buffet 
prediction accuracies necessary for use in final design. 

Basic research into the instability/breakdown mechanisms of concentrated vortices and free shear 
layers is needed 

- There is a lack of understanding of the flow physics for high angle, high subsonic speed 
conditions. Off-the-surface transonic features, such as terminating center-line shocks and shock- 
vortex structures, are known to occur and to have strong influences. 

Computationa l Vortex and Buffet Studies - Among a number of computational issues which could be 
highughted, one of the most .important is the establishment of required levels of algorithm accuracy 
and grid densities for nondissipative simulation of streamwise/leading-edge vortices Another issue 
which will be difficult to resolve involves the use of computational buffet predictions in the design 
process. How will such codes be validated? What will be the stated level of accuracy (or uncertainty?) 
m predicting buffeting response spectra and fatigue lifetimes? Experimental buffet data is usually 
available in the form of statistical measures. It will be interesting to note the sensitivity and 
repeatability of the corresponding computed statistical buffet measures to small changes in parameters. 


SUMMARY 

Computational methods and computer resource requirements for the prediction of tail buffeting have 
been surveyed. Computational finite difference methods have emerged which can provide accurate 
solutions of the thin-layer Navier-Stokes equations for flows about complete aircraft. Complex 
geometries have been modeled for low angles using 500,000 to 10 6 grid points for one-half span and 
requiring 5-40 hours of computer runtime per case. For higher angle, vortex-dominated flows 
computations for idealized geometries using 250,000 to 10^ grid points for one-half span require 2-25 
hours. The existing capability of CFD calculations of vortex-dominated flows offers encouragement 
with computed primary and secondary vortex locations and surface pressure levels in reasonably good 
agreement with experiment. Indications that vortex breakdown can be predicted are also encouraging 
although there have been no studies to verify the accuracy of computations of buffeting frequency 
response spectra or spatial correlations. All cases surveyed reported the need for further grid 
refinement in order to achieve quantitative agreement with experiment. Estimates of the number of grid 
points and the number the computational steps needed to perform buffet calculations have been 
developed. For current thin-layer Navier-Stokes codes using a full span grid of over 2 x 10^ grid 
points to model the F-15 airplane, approximately 1000 hours of computer runtime would be required 
for one buffet response calculation. Factors which could significantly lower or raise this estimate were 
discussed. Basic experimental studies are needed for the development of free shear layer turbulence 
models and for understanding of the fluid dynamics involved in the breakdown of concentrated 
vortices. Research is needed to develop numerical algorithms for convecting vortical flows which are 

nondissipative or which have controlled levels of dissipation such that accurate buffet calculations will 
be possible. 
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It can be anticipated that the three order of magnitude decrease in computer runtime needed to bring 
computational buffet predictions within reason will be forthcoming due to increases in computer speed 
and memory and increases in algorithm efficiency. A large unknown factor involves the nature of the 
origination and evolution of the buffet flowfield. It can also be anticipated that nonlinear fluid dynamic 
processes will be critical in the accurate simulation of such flows. This may lead to high sensitivity to 
initial conditions, with the result that estimates of runtimes and the required number of cases may be 
seriously underestimated. 
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Vortical Flow Regimes as a Function of Angle of Attack (from Reference 15) 
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Table n 


Computer Time Requirements for Buffet Calculations 


T cp „ = N„N„t 


X = 10 - 60 nsec /grid point /lime step 
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Table III - CFD Applications to Complex Aircraft Geometries 
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Table IV - CFD Calculations for Highly-Swept Wings at High Angles of Attack 
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- Euler code 

- sharp leading-edge 

- round leading-edge 

- vortex burst/unsteady calculations 
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Sharp leading-edge extension causes flow separation and vortex formation at high angles 



Figure 3. F-18 aircraft at high angle buffeting condition 
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b. Steady normal force coefficient versus angle-of-attack 
Figure 5. Vortex flow characteristics for a straked delta wing 6 






(a) u/Voo- 




a. contours of mean velocity 





(c) ui/Voo, percent. 

b. RMS velocity contours 


Figure 6. Three-component laser velocimetry flow survey for a YF-17 model in a low 
speed wind tunnel, a = 25 deg, Re = 326,000 J 1 
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Figure 7* Laser velocimetry flow survey of for an F-15 model in a low speed wind 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



Figure 8. Example of a complex separated-flow situation at deep buffet penetration 
conditions for a transonic speed. 



Figure 9. Development of a boundary layer on a flat plate. 
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SPIRAL TYPE VORTEX BREAKDOWN 



Figure 13. Schematic representation of spiral- type 
breakdown; longitudinal cross section.^ 


Region of strong fluctuations Zero flow line 



Figure 14. Diagram of vortex breakdown region JO 



a) Flood lamp 


b) Laser light sheet 


Figure 15. Smoke flow vizualization of leading-edge vortices on a shart-edged 85 deg. 
delta wing at a = 40 degJ^ 


37 



38 



(f) to” 40* cropped double* delta win a 


* 

t 



a. a = 12 deg. 


b. a = 20 deg. 


Figure 1 8. Top views of off-surface particle traces for three wing planforms.28 
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